Decorated Shastry-Sutherland lattice in the spin-^ magnet CdCu2(B03)2 
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We report the microscopic magnetic model for the spin-| Heisenberg system CdCu2(B03)2, one 
of the few quantum magnets showing the i-magnetization plateau. Recent neutron diffraction 
experiments on this compound [M. Hase et al, Phys. Rev. B 80, 104405 (2009)] evidenced long- 
range magnetic order, inconsistent with the previously suggested phenomenological magnetic model 
of isolated dimers and spin chains. Based on extensive density-functional theory band structure 
calculations, exact diagonalizations, quantum Monte Carlo simulations, third-order perturbation 
theory, as well as high-field magnetization measurements, we find that the magnetic properties of 
CdCu2(B03)2 are accounted for by a frustrated quasi-2D magnetic model featuring four inequivalent 
exchange couplings: the leading antiferromagnetic coupling Jd within the structural CU2O6 dimers, 
two interdimer couplings Jti and Jt2, forming magnetic tetramers, and a ferromagnetic coupling Jit 
between the tetramers. Based on comparison to the experimental data, we evaluate the ratios of 
the leading couplings Jd : Jti : Jt2 : Jit = 1 : 0.20 : 0.45 : —0.30, with Jd of about 178 K. The inequiva- 
lence of Jti and Jt2 largely lifts the frustration and triggers long-range antiferromagnetic ordering. 
The proposed model accounts correctly for the different magnetic moments localized on structurally 
inequivalent Cu atoms in the ground-state magnetic configuration. We extensively analyze the mag- 
netic properties of this model, including a detailed description of the magnetically ordered ground 
state and its evolution in magnetic field with particular emphasis on the ^-magnetization plateau. 
Our results establish remarkable analogies to the Shastry-Sutherland model of SrCu2(B03)2, and 
characterize the closely related CdCu2(B03)2 as a material realization for the spin-i decorated 
anisotropic Shastry-Sutherland lattice. 

PACS numbers: 71.20.Ps, 75.10.Jm, 75.40.Cx, 75.60.Ej 



I. INTRODUCTION 

Among the great variety of magnetic ground states 
(GSs), only few can be rigorously characterized using 
analytical considerations. To these exceptional cases 
belong, e.g., the Heisenberg zigzag chain model at the 
so-called Majumdar-Ghosh point, ^ a class of Kitaev 
models,^ as well as the dimerized phase of the Heisen- 
berg model on the Shastry-Sutherland lattice.'^ The lat- 
ter describes the magnetism of 5*= 1/2 spins on a square 
lattice, with two inequivalent exchange couplings: (i) 
a dimer-like coupling J connecting the spins pairwise 
along the diagonals of the square lattice, and (ii) the 
coupling J' along the edges of the square lattice (Fig. 1, 
left). For the ratio of antiferromagnetic (AFM) couplings 
J'/ J < 0.7, the GS is an exact product of singlets residing 
on the J bonds. ^ 

The Shastry-Sutherland model was first studied 
theoretically,'^ and only two decades later the spin S = ^ 
magnet SrCu2 (603)2 was claimed to be an experimen- 
tal realization of this model. ^ In the following years, nu- 
merous experimental and theoretical studies disclosed a 
rather complex magnetic behavior of this compound, in- 
cluding plateaux at 1/8, 1/4, and 1/3 of the satura- 
tion magnetization (Ref. 6) and an intricate coupling 
between magnetic and lattice degrees of freedom. For 
the experimentally defined ratio of the leading couplings 
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FIG. 1. (Color online) Left panel: the original Shastry- 
Sutherland model with the dimer-like coupling J (thick lines) 
and the coupling J' (thin lines) along the edges of the square 
lattice. Right panel: the decorated anisotropic Shastry- 
Sutherland model in CdCu2(B03)2. The horizontal and ver- 
tical couplings (J' in the original model) are not equivalent, 
leading to Jti (double line) and Jt2 (single line). In addition, 
the coupling Jti (dashed line) tiles the lattice into four-spin 
units (tetramers). Filled circles denote the decorating spins. 



J'/J = 0.64 (Refs. 8 and 9), the system is already in 
the dimerized singlet phase, but very close to the quan- 
tum critical point. ^'^ To access the critical point, further 
tuning of the system is required. For example, exter- 
nal pressure leads to at least one phase transition at 
low temperatures,^^ as evidenced by nuclear magnetic 
resonance (NMR) studies, while a recent theoretical 
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modcl^"^ with two inequivalcnt dimcr couplings J pre- 
dicts an unanticipated, effectively one-dimensional Hal- 
dane phase. ^'^'^^ 

An alternative way to tune the magnetic couplings is 
the directed substitution of nonmagnetic atoms or an- 
ionic groups. The impact of this substitution on the mag- 
netic coupling regime can be very different. The simplest 
effect is a bare change of the energy scale, which does not 
affect the physics, as for instance, in the spin-^ uniform- 
chain compounds Sr2Cu(P04)2 and Ba2Cu(P04)2.^^ An- 
other possibility is alteration of the leading couplings, as 
in the frustrated square lattice systems AA'YO(P04)2 
(Ref. 16) or the kagome lattice compounds kapellasite 
Cu3Zn(OH)6Cl2 and haydeeite Cu3Mg(OH)6Cl2.^^ Such 
systems have excellent potential for exploring magnetic 
phase diagrams. Finally, in some cases, chemical substi- 
tutions drastically alter the underlying magnetic model, 
as shown recently for the family of CU2A2O7 compounds 
(A = P,As,V).i8 

However, the most common effect of chemical substi- 
tution is a structural transformation, which apparently 
leads to a major qualitative change in the magnetic cou- 
plings. For instance, the Se-Te substitution, accompa- 
nied by a major structural reorganization, transforms the 
Heisenberg-chain system CuSe205 (Ref. 19) into the spin 
dimer compound CuTe205 (Ref. 20). Likewise, the La- 
Bi substitution in the 2D square lattice system La2Cu04 
(Ref. 21) leads to the three-dimensional magnet Bi2Cu04 
(Ref. 22), etc. 

Coming back to the S = Shastry-Sutherland com- 
pound SrCu2(B03)2, the directed substitution of non- 
magnetic structural elements is an appealing tool to ex- 
plore the phase space of the Shastry-Sutherland model. ^'^ 
One possible way is the replacement of Sr^"'' by another 
divalent cation. Indeed, the compound CdCu2(B03)2 
exists, and its magnetic properties have been investi- 
gated experimentally.^^"^^ In contrast to SrCu2(B03)2, 
the magnetic GS of CdCu2(B03)2 is long-range ordered, 
while the magnetization M{H) curve exhibits a plateau 
at one- half of the saturation magnetization.^^ 

The substitution of Sr by Cd in SrCu2(B03)2 is ac- 
companied by a structural transformation. Thus, the 
initial tetragonal symmetry is reduced to monoclinic, en- 
gendering two independent magnetic sites: Cu(l) and 
Cu(2). The presence of Cu(l) residing in the structural 
dimers, similar to SrCu2(B03)2, and the relatively short 
Cu(2)-Cu(2) connections along c motivated the authors 
of Ref. 24 to advance a tentative magnetic model com- 
prising isolated magnetic dimers of Cu(l) pairs and in- 
finite chains built by Cu(2) spins. This model was in 
agreement with the experimental magnetic susceptibil- 
ity x(2^)j magnetization M{H), and electron spin reso- 
nance (ESR) data. However, later neutron diffraction 
(ND) studies challenged this model, since the ordered 
magnetic moments localized on Cu(l) amount to 0.45 /xb 
compared to 0.83 /xb on Cu(2).26 Therefore, the Cu(l) 
dimers do not form a singlet state with zero ordered mo- 
ment, but rather experience a sizable coupling to Cu(2). 



This coupling is ferromagnetic, as inferred from the mag- 
netic structure, and contrasts with the AFM coupling 
assumed in the previous studies. 

The tentative magnetic model from Ref. 24 was 
founded on the assumption that the strongest magnetic 
couplings are associated with the shortest Cu-Cu inter- 
atomic distances. However, such an approach fully ne- 
glects the mutual orientation of magnetic orbitals, which 
play a decisive role for the magnetic coupling regime. 
Moreover, the experimental results suggest that the ac- 
tual spin model of CdCu2(B03)2 is more complex than 
the simple "chain -I- dimer" scenario. 

In this paper, we investigate the magnetic prop- 
erties of CdCu2(B03)2 on a microscopic level, using 
density-functional-theory (DFT) band structure calcula- 
tions, exact diagonalization (ED) and quantum Monte 
Carlo (QMC) numerical simulations, as well as ana- 
lytical low-energy perturbative expansions and linear 
spin-wave theory. We show that the magnetism of 
CdCu2(B03)2 can be consistently described by a two- 
dimensional (2D) frustrated isotropic model (crystalline 
and exchange anisotropy effects are neglected) with four 
inequivalcnt exchange couplings, topologically equivalent 
to the decorated anisotropic Shastry-Sutherland lattice 
(see Fig. 1, right). Therefore, the chemical similarity 
between SrCu2(B03)2 and CdCu2(B03)2 is retained on 
the microscopic level. The dominance of the exchange 
interaction within the structural Cu(l) dimers al- 
lows to describe the magnetic properties within an ef- 
fective low-energy model for the Cu(2) sites only. This 
effective model explains the sizable staggered magneti- 
zation of Cu(l) dimers despite the large Jd, and unrav- 
els the nature of the wide 1/2-magnetization plateau in 
CdCu2(B03)2. The experimental data, both original and 
taken from the literature, are reconsidered in terms of the 
DFT-based microscopic model and the effective model. 
As a result, we find an excellent agreement between the 
simulated and measured quantities. 

This paper is organized as follows. Sec. II is de- 
voted to the discussion of the structural peculiarities 
of CdCu2(B03)2. The methods used in this study are 
presented in Sec. III. Then, the microscopic magnetic 
model is evaluated using DFT band structure calcula- 
tions, and its parameters are refined using ED fits to 
the experiments (Sec. IV). The magnetic properties of 
CdCu2(B03)2 are in the focus of Sec. V. We show that 
the magnetism of CdCu2(B03)2 is captured by an ef- 
fective low-energy description of the Cu(2) spins (on the 
Shastry-Sutherland lattice) , which are further studied us- 
ing linear spin- wave theory as well as QMC, and corrobo- 
rated by ED for the microscopic magnetic model. Finally, 
in Sec. VI a summary and a short outlook are given. 



II. CRYSTAL STRUCTURE 

The low-temperature crystal structure of 
CdCu2(B03)2 is reliably established based on syn- 
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FIG. 2. (Color online) The crystal structure of CdCu2(B03)2 
consists of structural Cu(l)206 dimers, Cu(2)04 distorted 
plaquettes, BO3 triangles connected pairwise as well as iso- 
lated Cd atoms (spheres). Oxygen atoms are shown as small 
spheres located at the vertices of Cu(l)04, Cu(2)04 and BO3 
polyhedra. The numbers "1" and "2" denote Cu(l) and 
Cu(2), respectively. A magnetic tetramer (see Sec. IV A) is 
shaped by a gray oval. Right bottom: position of the magnetic 
layers (shaded) in the crystal structure of CdCu2(B03)2. 



chrotron x-ray and neutron diffraction data.^^ The basic 
structural elements of CdCu2 (603)2 are magnetic layers 
(Fig. 2) that stretch almost parallel to (102) planes 
(Fig. 2, right bottom). These layers are formed by Cu-0 
polyhedra coupled by B2O5 units. 

Although Cd^"'' is isovalent to Sr^"*", the substan- 
tial difference in their ionic radii (1.26 A and 0.95 A, 
respectively)^'' gives rise to different structural mo- 
tives in the two systems. Thus, in sharp contrast 
with SrCu2(B03)2, the monoclinic crystal structure of 
CdCu2 (603)2 comprises two inequivalent positions for 
Cu atoms: Cu(l) and Cu(2). The local environment of 
these two sites is different: Cu(l) forms Cu(l)04 plaque- 
ttes that in turn build structural Cu(l)206 dimers, re- 
sembling SrCu2 (603)2, while Cu(2) atoms form tctrahc- 
drally distorted Cu(2)04 plaquettes that share a common 
O atom with the Cu(l)206 dimers (Fig. 2). As a result, 
structural Cu(l)2Cu(2)20i2 tetramers are formed. In- 
terestingly, the structural tetramers do not coincide with 
the magnetic tetramers (shaded oval in Fig. 2) that are 
the basic units of the spin lattice (Sec. IV A). 

Large cavities between the neighboring magnetic 
tetramers accommodate Cd atoms. As a result, the space 
between the magnetic layers is rather narrow, leading 
to short interatomic distances (^3.4 A) between the Cu 
atoms belonging to neighboring magnetic layers. On the 
basis of such short Cu-Cu distance, the relevance of the 



respective interlaycr coupling (the chain-like coupling J3 
in the notation of Ref. 24) has been conjectured.^'* How- 
ever, the orientation of magnetic Cu(l)206 and Cu(2)04 
polyhedra (almost coplanar to the magnetic layers) ex- 
cludes an appreciable interlayer coupling. This quasi-2D 
scenario is readily confirmed by our DFT calculations 
(Sec. IV A). 

The evaluation of the intralayer coupling based on 
structural considerations is more involved. In many 
cupratcs, the leading exchange couplings are typically 
associated with the Cu-Opiaq-Cu connections, i.e., when 
the neighbouring plaquettes share a common corner (one 
Opiaq atom) or a common edge (two Opiaq atoms). ^® In 
the crystal structure of CdCu2(603)2, both types of con- 
nections arc present. 

To estimate the respective magnetic exchange cou- 
plings, the Goodenough-Kanamori rules^^ are typically 
applied. Using a closely related approach, 6raden et a/.^° 
evaluated the dependence of the magnetic exchange J on 
the Cu-O-Cu angle. According to their results, for the 
Cu~0-Cu angle of 90°, the resulting total coupling J, 
being a sum of AFM J^™ (positive) and FM J^^ (neg- 
ative) contributions, is FM. An increase in this angle (i) 
decreases the contribution caused by Hund's cou- 
pling on the O sites, ■^^ and (ii) enhances the Cu-O-Cu 
superexchange which in turn increases 

JAFM^ As a re- 

suh, for the Cu-O-Cu angle close to 97°, J^^ and J"^^^ 
are balanced, thereby cancelling the total magnetic ex- 
change (J = 0). Further increase in the Cu-O-Cu angle 
breaks this subtle balance and gives rise to an overall 
AFM exchange, which grows monotonically up to Cu- 
O-Cu =180°. 

Here, we try to apply this phenomenological approach 
and estimate the respective exchange integrals. The in- 
tradimcr Cu(l)-0-Cu(l) angle in CdCu2 (603)2 is 98.2°, 
only slightly exceeding 97°. Therefore, a weak AFM ex- 
change coupling could be expected. For the Cu(l)-0- 
Cu(2) corner-sharing connections, the angle amounts to 
117.5°, which should give rise to a sizable AFM exchange. 
However, our DFT+U calculations deliver a magnetic 
model, which is in sharp contrast with these expecta- 
tions: the AFM coupling within the structural dimers 
is sizable, while the coupling along the Cu(l)-0-Cu(2) 
connections is much smaller and FM. The substantial dis- 
crepancy between the simplified picture, based on struc- 
tural considerations only, and the microscopic model, evi- 
dences the crucial importance of a microscopic insight for 
structurally intricate magnets like CdCu2 (603)2. 



III. METHODS 

DFT band structure calculations were performed us- 
ing the full-potential code fplo-8. 50-32.^^ For the struc- 
tural input, we used the atomic positions based on the 
neutron diffraction data at 1.5 K: space group P2i/c (14), 
a = 3.4047A, &=15.1376A, c = 9.2958A, /3 = 92.807°. 
For the scalar-relativistic calculations, the local-density 
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approximation (LDA) parameterization of Perdew and 
Wang has been chosen. '^'^ We also cross-checked our re- 
sults by using the parameterization of Perdew, Burke 
and Ernzerhof based on the generalized gradient approx- 
imation (GGA)."^^ For the nonmagnetic calculations, a 
/cnicsh of 30x6x12 = 2160 points (728 points in the ir- 
reducible wedge) has been adopted. Wannier functions 
(WF) for the Cu 3dx2-y2 states were evaluated using 
the procedure described in Ref. 35. For the supercell 
magnetic DFT-|-?7 calculations, two types of supercells 
were used: a supercell metrically equivalent to the unit 
cell (sp. gr. PI, 4x2x2 ftmesh) and a supercell doubled 
along c (sp. gr. PI, 4x1x1 fcmesh). Depending on the 
double-counting correction (DCC),''^ the around- mean- 
field (AMF) or the fully localized limit (FLL), we var- 
ied the on-site Coulomb repulsion parameter U^d in a- 
wide range: Uad = 5.5-7.5 eV for the AMF and U^d = 8.5- 
10.5 eV for the FLL calculations,'^^ keeping J3£j = leV. 
The local spin-density approximation (LSDA)-I-C/ results 
were cross-checked by GGA+U calculations. All results 
were accurately checked for convergence with respect to 
the /cmesh. 

Full diagonalizations of the resulting Heisenberg 
Hamiltonian were performed using the software pack- 
age ALPS-1.3 (Ref. 38) on the N=16 finite lattice with 
periodic boundary conditions. Translational symmetries 
have been used (the code does not allow for rotational 
symmetries). Spin correlations in the GS as well as the 
lowest-lying excitations were computed by Lanczos diag- 
onalization of = 32 sites finite lattices using the code 
SPINPACK.'^^ Translational symmetries and the commu- 
tation relation [H, Sz] = were used. The finite lattices 
are depicted in Fig. 12. The QMC simulations of the ef- 
fective model were performed on finite lattices with up 
to TV = 2304 spins with periodic boundary conditions us- 
ing the looper*^ algorithm from the ALPS package. '^'^ For 
T = 0.05 J ~ 0.0031 Jd(~ 0.55 K), we used 40 000 loops for 
thermalization and 400 000 loops after thermalization. 

A powder sample of CdCu2 (603)2 was prepared by an- 
nealing a stoichiometric mixture of CuO, CdO, and B2O3 
in air at 750 °C for 96 hours. The sample contained trace 
amounts of unreacted CuO and Cd2B205 (below 1 wt.% 
according to the Rietveld refinement), as evidenced by 
powder x-ray diffraction (Huber G670 camera, CuKq,i ra- 
diation, ImagePlate detector, 26* = 3-100° angle range). 
To improve the sample quality, we performed additional 
annealings that, however, resulted in a partial decompo- 
sition of the target CdCu2(B03)2 phase. Both impurity 
phases reveal a low and nearly temperature-independent 
magnetization within the temperature range under inves- 
tigation. Therefore, they do not affect any of the results 
presented below. 

The magnetic susceptibility (x) was measured with an 
MPMS SQUID magnetometer in the temperature range 
2 — 380 K in applied fields up to 5 T. The high-field mag- 
netization curve was measured in pulsed magnetic fields 
up to 60 T at a constant temperature of 1.5 K using the 
magnetometer installed at the Dresden High Magnetic 



Field Laboratory. Details of the experimental procedure 
can be found in Ref. 41. 



IV. EVALUATION OF A MICROSCOPIC 
MAGNETIC MODEL 

A. Band structure calculations 

We begin our analysis with the evaluation of the mag- 
netically relevant orbitals and couplings. The LDA yields 
a valence band with a width of about 9 eV, dominated 
by O 2p states (Fig. 3, top panel). Interestingly, the 
lower part of the valence band exhibits sizable contri- 
bution from the Cd 4d states, hinting at an appreciable 
covalcncy of Cd-0 bonds. The energy range between 
—3 and 0.5 eV is dominated by Cu 3d states. The pres- 
ence of bands crossing the Fermi level evidences a metal- 
lic GS, which contrasts with the emerald-green color of 
CdCu2(B03)2. This well-known shortcoming of the LDA 
and the GGA originates from underestimation of strong 
correlations, intrinsic for the electronic configuration 
of Cu^"''. To remedy this drawback, we use two ap- 
proaches accounting for the missing part of correlations: 
(i) on the model level, via mapping the LDA bands onto 
a Hubbard model (model approach); and (ii) directly 
within a DFT code, by adding an energy penalty Usd 
for two electrons occupying the same orbital (DFT-|-?7 
approach) . 

energy (eV) 
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FIG. 3. (Color online) Top: LDA valence band of 
CdCu2(B03)2. Atom-resolved DOS (number of states per eV 
and crystallographic unit cell) are shown. Bottom left: LDA 
band structure of the well-separated eight-band complex at 
the Fermi level and the WF tight-binding fit within an effec- 
tive one-orbital model. Bottom right: atom-resolved DOS of 
the eight-band complex. 

The states relevant for magnetism are located in a close 
vicinity of the Fermi level. Common for cuprates, these 
states form a band complex, well-separated from the rest 



of the valence band. This band complex, depicted in 
Fig. 3 (bottom panels), comprises eight bands, in accord 
with the eight Cu atoms in a unit cell. The projection 
of these states onto a local coordinate system^^ readily 
yields their orbital character: typical for cuprates, these 
states correspond to the antibonding cr-combination of 
Cu 3^3.2 _j,2 and O 2p orbitals. To visualize this combi- 
nation in real space, the WFs for Cu 3dx2-y2 states can 
be used (Fig. 4). 

To evaluate the relevant couplings, we calculate hop- 
ping matrix elements of the WFs and introduce them as 
parameters {U) into an effective one-orbital tight-binding 
model. This parameterization is justified by an excellent 
fit (Fig. 3, bottom left) to the LDA band dispersions. The 
relevant couplings (t^ > 20 meV) are presented in Table I 
(third column). 

The transfer integrals ti can be subdivided into four 
groups, according to their strength. The first group con- 
tains only one term td ("d" stands for "dimer"), which 
clearly dominates over all other couplings, despite the 
respective Cu(l)-0-Cu(l) angle being rather close to 
90°. The second group comprises two terms, tti and 
tt2 ("t" stands for "tetramer") that couple Cu(l) and 
Cu(2) atoms via Cu(l)-0-0-Cu(2) paths (Fig. 4) within 
the magnetic tetramers. Although the Cu(l)-Cu(2) dis- 
tance for tti is shorter, the respective coupling is slightly 
smaller than it2- This difference plays an important role 
for the magnetism of CdCu2 (603)2, as will be shown 
later. 

The couplings t^, tji, and tt2 form magnetic tetramers 
(shaded oval in Fig. 2). These tetramers are coupled to 
each other by two interactions forming the third group: 
tit and tcd- The coupling ^it ('it" stands for "intcrte- 
tramer") corresponds to the magnetic exchange J2 be- 
tween dimers and chains in the notation of Ref. 24. It 
is short-ranged, and runs via the corner-sharing Cu(l)- 
0-Cu(2) connection. The second intertetramer term tcd 
couples Cu(l) and Cu(2) via Cd atoms (not shown). Fi- 
nally, the fourth group of transfer integrals includes all 
other terms that are smaller than 20 meV. Since these lat- 
ter terms play a minor role for the magnetic properties 
of CdCu2(B03)2, they are disregarded in the following 
discussion. 

To account for electronic correlations, neglected in the 
tight-biding model, the transfer integrals ti are mapped 
onto a Hubbard model. This way, under the condi- 
tions of half-filling and strong correlations (UeS^ti), 
the AFM contribution J^^^ to the exchange integrals 
can be estimated in second-order perturbation theory:^'^ 
jAFM _^^2^^^^_44 jjqpq^ jg g^jj effective term which 

takes into account the correlations in the magnetic Cu 
5dj.2_y2 orbitals, screened by the O 2p orbitals (their 
strong hybridization can be seen in Fig. 4). Although 
the experimental evaluation of UeS is challenging (see, 
e.g., Ref. 45), there is an empirical evidence that UcS 
values in the 4-5 eV range yield reasonable results for 
Cu^"*" compounds. ^^'^^'^^ Here, we adopt C/off = 4.5eV to 
estimate jf^^^^ (Table I, fourth column). 
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FIG. 4. (Color online) Wannier functions for the Cu 3d^2_y2 
states. Light and dark areas (colors) denote different phases 
of the functions. Relevant exchange couplings are depicted 
by lines: thick (red-black) lines (Jd), thinner solid (Jt2) and 
double (Jti) lines, as well as dashed lines (Jit). Note the 
Cu-O-O-Cu superexchange paths via the BO3 triangles. 



TABLE I. Cu-Cu distances (rf, in A), transfer integrals {ti, 
in meV), antiferromagnetic jf^^^ (in K) contributions, as well 
as the total exchange integrals Ji (in K) for the leading cou- 
plings in CdCu2 (303)2. The values of Ji are evaluated within 
the DFT+U method for two different double- counting correc- 
tion schemes with Usd = 6.5 eV for AMF and 9.5 eV for FLL. 



exchange 


d 


t^ 


jAFM 


J, (LSDA+(7) 
AMF FLL 


Jd 


2.957 


203 


425 


146 


202 


Jti 


5.268 


73 


55 


30 


29 


Jt2 


6.436 


85 


75 


64 


55 


Jit 


3.228 


47 


23 


-45 


-135 


Jcd 


6.450 


41 


17 


5 


6 


Jti : Jd (%) 
Jit : Jd (%) 
Jti : Jt2 (%) 






20 
-31 
47 


14 
-67 
53 



Resorting from ti to Ji basically preserves the pic- 
ture of coupled tetramers with a strong dimer-like cou- 
pling. However, the interatomic Cu-Cu distances for Jd 
and Jit are rather small, thus sizable FM contributions 
to these couplings can be expected. To account for the 
FM exchange, we perform DFT+U calculations for dif- 
ferent magnetic superccUs,^^ and subsequently map the 
resulting total energies onto a Heisenberg model. This 
way, the FM contribution can be evaluated by subtract- 
ing j/^^^ estimates from the values of the total exchange 
Ji, obtained from DFT+U calculations. The results of 
DFT+U calculations typically depend on (i) the func- 
tional (LSDA or CCA), (ii) the DCC scheme (AMF or 
FLL). and (iii) the value of Usd used.^^ Representative 
values for the total exchange integrals Ji arc summarized 
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in Tabic I. 

We first discuss common trends that do not depend on 
a particular DFT+U calculational scheme. First, despite 
the substantial FM contribution, Jd preserves its role of 
the leading coupling. Second, the FM contributions fur- 
ther enhance the difference between Jti and Jt2- This 
difference is crucial, since these two couplings are frus- 
trated, and the degree of frustration largely depends on 
their ratio. Third, the interdimcr coupling Jit is clearly 
FM, while a sizable FM contribution to Jcd practically 
switches this coupling off. 

The DFT+U results exhibit a rather strong depen- 
dence on the calculational scheme and the parameters 
used. A comparative analysis of the exchange values 
shows that the couplings Jd, Jti, Jt2, and Jcd yielded 
by GGA-I-J7 are reduced by 20-30% compared to the 
LSDA-ft/ estimates (for the same DCC and C/sd). The 
difference between the LSDA-I-J7 and GGA+U results for 
Jit is slightly larger. However, the ratios of the leading 
couplings, which determine the nature of the GS, are very 
similar for both functionals. 

The DCC scheme has a much larger impact on the 
values of exchange integrals. In particular, FLL yields 
substantially larger Ja values than AMF. In addition, 
the Jti/ Jd and Jt2/Jd ratios are slightly reduced in FLL 
compared to the AMF estimates (Table I). 

The dependence of the DFT-|-?7 results on Usd shows 
a clear and robust trend: larger Usd lead to smaller abso- 
lute values of the leading exchange couplings. However, 
the slope of the Ji{Uzd) dependence (not shown) slightly 
differs for different couplings. As a result, the ratios of 
the couplings also depend moderately on the C/a^ value. 

In conclusion, DFT provides a robust qualitative mi- 
croscopic magnetic model and unambiguously yields the 
signs of the leading Jj 's, as revealed by excellent agree- 
ment with the experimentally observed magnetic struc- 
ture (Fig. 5). However, the accuracy of the individual 
magnetic couplings is not enough to provide a quantita- 
tive agreement with the experiments. In the next step, 
the evaluated model is refined by a simulation of its mag- 
netic properties and subsequent comparison to the exper- 
iments. 



B. Simulations 

Since the microscopic magnetic model is 2D and 
frustrated, numerically efficient techniques that account 
for the thermodynamic limit, such as quantum Monte 
Carlo (QMC) or density-matrix renormalization group 
(DMRG), can not be applied. Although new meth- 
ods are presently under active development,^* a stan- 
dard method to address such problems is exact diag- 
onalization (ED) of the respective Hamiltonian matrix 
performed on finite lattices of N spins. '^^ However, the 
performance of present-day computational facilities lim- 
its the size of feasible finite lattices to = 24 sites^° for 
finite-temperature properties or A^ = 42 sites^^ for the 




FIG. 5. (Color online) Microscopic magnetic model for 
CdCu2(B03)2. The leading magnetic couplings are denoted 
by lines. The experimental magnetic structure (Ref. 26) is 
depicted by arrows. The crystal structure sketch is the same 
as in Fig. 2. 



GS. Additional restrictions arise from topological fea- 
tures of the particular model, e.g., only certain values 
of N fit to periodic boundary conditions. Here, we apply 
ED to study the magnetic properties of CdCu2 (603)2 
using exact diagonalization on N =1Q sites (for the mag- 
netic susceptibility) and A^ = 32 sites (for GS spin-spin 
correlations and the GS magnetization process in mag- 
netic field). 

Due to the ambiguous choice of the DCC and U^id, 
our DFT calculations do not provide a precise position 
of CdCu2 (603)2 in the parameter space of the proposed 
Jd^Jti'^Jt2^Jit model. However, the ratios of the leading 
couplings (Table I, three bottom lines) follow a distinct 
trend: the AMF solutions yield larger Jti : Jd and sub- 
stantially smaller |Jit|: Jd values, than FLL. Compared 
to the difference between the AMF and FLL solutions, 
the value of U^d and the functional used influence the 
ratios only slightly. Therefore, we adopt the two differ- 
ent DFT-t-C/ solutions: the AMF solution Jd : Jti : Jt2 : -'^it 
~ 1 : 0.20 : 0.45 : -0.30 (AMF, = 6.5 eV) and the FLL 
solution Jd : Jti : Jt2 : Jit 1 : 0.15 : 0.25 : -0.70 (FLL, 
t^3d = 9.5eV) as starting points for simulations of mag- 
netic susceptibility. 

We have simulated the magnetic susceptibility on a 
A^= 16 sites finite lattice (Fig. 12, bottom) with periodic 
boundary conditions, and by keeping the Jd : Jti : Jt2 '■ Jit 
ratios fixed. Such simulations yield the reduced mag- 
netic susceptibility x*(T*), which can be compared to 
the experimental x(2^) by fitting free parameters: the 
overall energy scale Jd, the g-factor and the temperature- 
independent contribution Xo-^^ In accord with the 
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FIG. 6. (Color online) Top panel: ED (iV = 16 sites) fit to 
the experimental (exp.) magnetic susceptibility. Deviations 
at r<30K are an extrinsic finite-size effect. Inset: fit to 
the low-temperature xiT) using QMC for the effective model 
(Sec. VA) on a N =24x24 sites finite lattice. Note the im- 
proved agreement between 15 and 30 K, compared to the ED 
fit. Bottom panel: magnetization curve measured at 1.5 K 
and simulated by ED (finite-temperature magnetization for 
A'' = 16 and the GS magnetization for N = 32 sites) scaled us- 
ing the parameters (Jd and g-factor) obtained from the fit to 
X{T) (Table II, top line). 



Mermin- Wagner theorem, ^'^ our 2D magnetic model does 
not account for long-range magnetic ordering, hence 
the simulated magnetic susceptibility is valid only for 
T>Tn = 9.8K.^^ Moreover, at low temperatures (even 
above Tn) finite-size effects become relevant. To find 
an appropriate temperature range, we used an auxil- 
iary eight-order high-temperature series expansion, and 
thus evaluated 30 K as a reasonable lowest fitting temper- 
ature for the N= 16 sites finite lattice. 

The fitted values of Jd, g, and xo ^-re summarized in 
Table II. Since the FLL solution (Table II, bottom line) 
yields unrealistically large g = 2.457, it can be safely ruled 
out. By contrast, the AMF solution (Table II, first row) 
yields a reasonable value of y = 2.175.^* The resulting fit 
to the experimental xiT) curve is shown in Fig. 6 (top). 

Next, we challenge the solutions from Table II using 
the high-field magnetization data (Fig. 6, bottom). In 
particular, the critical field Hd, corresponding to the 
onset of the magnetization plateau (Fig. 6, bottom), can 
be estimated independently for both experimental and 
simulated curves. Experimentally, Hd = 23.2 T was eval- 
uated by a pronounced minimum in d'^m/dH'^.^^ The 
simulated magnetization curves m(h*)^ where h* is the 
reduced magnetic field in the units of Jd, were scaled to 
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FIG. 7. (Color online) The three topologically different clus- 
ters that contribute up to third order in the effective exchange 

Jij in the limit Ji,2,3 ^ Jd, where Ji, J2, and J3 belong to the 

(2) (3) 

set: Jtl, Jt2, or Jit. Here J^' and denote, respectively, 
the contributions from second and third-order perturbation 
theory. For any given cluster, we include only processes that 
involve each of the weak couplings at least once. 



the experimental data using the expression 

H = {kBJdg-'ti^^)h*, (1) 

and by adopting the values of Jd and g from the fits to 
x(T). This way, Hd can be estimated directly, i.e., with- 
out using additional parameters beyond Jd and g. Such 
Hd estimates are given in the last column of Table II. 
Both the AMF and the FLL solutions are in excellent 
agreement with the experimental Hd- Therefore, only 
the AMF solution Jd : Jti : Jt2 -Jit^l- 0.20 : 0.45 : -0.30 
yields good agreement for both the g-factor and the Hd 
values. Thus, in the following, we restrict ourselves to 
the analysis of this solution (Table II, first row). 

So far, the simulations evidence a good agreement be- 
tween the DFT-based model and the experimental in- 
formation on thermodynamical properties. In the fol- 
lowing, we perform an elaborate theoretical study of the 
proposed model. A deep insight into its physics enables 
further experimental verification of our microscopic sce- 
nario and allows to get prospects for new experiments. 

V. PROPERTIES OF THE MAGNETIC MODEL 

A. Effective low-energy theory 

From the above discussion, it has become clear that 
the intradimer exchange Jd is considerably larger than 



TABLE II. Fitted values of Jd (in K), the g-factor and xo (in 
10"'' X emu [mol Cu]~'), as well as the critical field Hd (in 
T), corresponding to the onset of a magnetization plateau. 
The simulations have been performed on A" = 16 sites finite 
lattices using ED. 
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FIG. 8. (Color online) The effective interactions between the Cu(2) sites (filled circles) up to lowest second order. Left panel; 
Eight paths a — h that provide a nonzero exchange interaction up to second order. Middle panel: effective model based on 
a — h paths: Jx = c -\- d, Jy = a + b, and J = e + f + g + h. In third order, one should also include the processes from the third 
cluster of Fig. 7, apart from a renormalization of the amplitudes a — h (see text). Right panel: The resulting effective model in 
the topologically equivalent square lattice version of the structure. The indices 1-4 label the four sites of the unit cell of this 
lattice model. The arrows denote the classical GS of the model. 



the remaining energy scales in the problem. This separa- 
tion of energy scales leads naturally to the idea that the 
physics of the problem can be understood by a pertur- 
bative expansion around the limit Jti = Jt2 = Jit = 0. In 
this limit, the Cu(l) sites coupled by Jd form quantum- 
mechanical singlets, while the remaining, Cu(2) sites arc 
free to point up or down, which defines a highly degener- 
ate GS manifold. By turning on the remaining couplings 
Jti, Jt2, and Jit, the Cu(2) sites begin to interact with 
each other through the virtual excitations of the Cu(l) 
dimers out of their singlet GS manifold. By integrating 
out these fluctuations, one can derive an effective low- 
energy model for the Cu(2) sites only. As we show be- 
low, a degenerate perturbation theory up to third order 
gives a remarkably accurate description of the low-energy 
physics of the problem. 

Figure 7 shows the three topologically different clus- 
ters G that involve a given Cu(l) dimcr and contribute 
up to third order in perturbation theory. To avoid dou- 
ble counting of the processes that live on subclusters of 
G,^^ wc keep only those processes that invoke each weak 
coupling of G at least once. Here the weak couplings 
■^1,2,3 <CJd belong to the set: Jti, Jt2, and Jit. 

For the following discussion, it is essential to empha- 
size the difference between the first two clusters in Fig. 7: 
in (a) the two couplings involve the same site a of the 
dimcr, while in (b) they involve both sites a and /?. Since 
the singlet wave function \s)ai3 = ^ iltaifs) - liatfs)) 
on the strong dimcr is antisymmetric with respect to in- 
terchanging a and /S, the two processes have opposite 
sign (and are equal in magnitude up to second order). In 
particular, if Ji J2 > 0, the effective coupling between the 
sites i and j is ferromagnetic in (a) but antiferromagnetic 
in (b). For the third cluster, we have a third-order pro- 
cess where each of the three weak couplings Ji , J2 and 
J3 has been involved once in the whole process. 



Using these quite general results, we may turn to our 
spin lattice and examine the possible effective processes 
that are active up to third order. We first discuss the 
contributions of the type shown in Fig. 7(a) and (b). 
Altogether, we find eight different paths a~h [these are 
shown in Fig. 8(a)] with the following amplitudes: 
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In the right-hand side equalities, we made use of the esti- 
mated ratios of the original couplings from Table II (the 
AMF solution). Including now the contributions of the 
type shown in Fig. 7(c), we get altogether three different 
effective exchange interactions between the Cu(2) spins: 



Jj. = c + d — 



Jti Jt2 Jit 



-0.05578Jd , 



(2) 
(3) 



Jy = a + h —2 — ~ -0.03234Jd , 

2Jj 

7 I f I I ; Jti Jt2( Jti + Jt2) nnci-TOT t a\ 

J^e + f + g + h —2 ~ 0.06172 Jd (4) 
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These couplings are shown in the right panel of Fig. 8, 
where we use the square-lattice representation that is 
topologically equivalent to our effective model. Hence, 
the resulting effective model is the anisotropic Shastry- 
Sutherland lattice model with an AFM exchange cou- 
pling J for the diagonal bonds, an AFM exchange Jx for 
the horizontal bonds, and an FM exchange Jy for the 
vertical bonds. This model readily explains the observed 
magnetic ordering on the Cu(2) sites since, on the classi- 
cal level, one can minimize all interactions simultaneously 
in the collinear "stripe" configuration shown in the right 
panel of Fig. 8. 

To challenge the applicability of the effective model, we 
refer again to the experimental x{T) dependence. The 
nonfrustrated nature of the effective model enables us to 
use numerically efficient QMC techniques.'^* To this end, 
we have simulated x*{T*) of the effective model on a fi- 
nite iV = 24x24 sites lattice with periodic boundary con- 
ditions. A parameter-free fit to the experimental curve 
is obtained by using Eqs. (2)-(4) above, and the values 
of Jd and g from the ED fit (Table II). By construction, 
the effective model is valid at low temperatures (T^ Jd), 
but above the long-range magnetic ordering temperature 
T/v. In this temperature range, the effective model yields 
excellent agreement with the experimental x(^) (see in- 
set in top panel of Fig. 6), justifying it as an appropriate 
low-energy model for CdCu2 (603)2. 



B. The size of the Cu(2) moments: linear spin wave 
theory in the effective model and QMC simulations 

The size of the Cu(2) moment in the GS of the system 
is not equal to the classical value gfis/^, but it will be 
somewhat lower due to quantum fluctuations. A similar 
correction is expected, e.g., for the GS energy. We can 
calculate these corrections by a standard 1/S semiclassi- 
cal expansion around the stripe phase keeping only the 
quadratic portion of the fluctuations. The details of this 
calculation are provided in Appendix A, and here we only 
discuss the main results. 

Figure 9 shows the magnon dispersions of the effec- 
tive model along certain symmetry directions in the Bril- 
louin zone (BZ). There are four magnon branches but 
these pair up into two 2-fold degenerate branches due to 
the symmetry of interchanging the sites (1,2) with (3,4) 
in the unit cell (see right panel of Fig. 8). The accu- 
racy of the magnon dispersions was cross-checked by the 
McPhase code.^*" 

Let us now look at the quadratic correction to the 
GS energy which, as shown in Appendix A. consists of 
two terms, 6E1 and (5i?2- The first term is given by 
SEi = -2^NucS, where ^ = J + 2Jx - 2Jy, and Nuc 
is the number of unit cells. So adding 5Ei to the classi- 
cal GS energy E'ciass = —'^S.NucS'^ , amount s to a renor- 
malization of the spin from 5* to S{S + 1). The second 
quadratic correction SE2 stands for the zero-point energy 
of the final, decoupled harmonic oscillators of the theory 
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FIG. 9. (Color online) Magnon dispersions (along certain 
symmetry directions in the BZ) obtained from linear spin 
wave theory around the stripe phase of the Cu(2) spins (see 
right panel of Fig. 8). There are four magnon branches but 
these pair up into two 2-fold degenerate branches, which is due 
to the symmetry of interchanging the sites (1,2) with (3,4) of 
the unit cell (see right panel of Fig. 8). 

[see Eq. (A13)], and can be calculated by a numerical 
integration over the BZ. Altogether, we flnd 

Eci,ss/Nuc = -0.119 Jd (5) 
SEi/Nuc = -0.238 Jd (6) 
(5S2/7V„c^ +0.2092 Jd. (7) 

Thus, quadratic fluctuations actually reduce the GS en- 
ergy by about 24%. We now turn to the renormaliza- 
tion of the magnetic moment, 5S. As explained in Ap- 
pendix A, this can also be calculated by a numerical in- 
tegration over the BZ, see Eq. (A16). We find 

(JS* -0.159^5' = 0.5 -^S*- 0.341 . (8) 

With g ~ 2.175 we get a local magnetic moment on the 
Cu(2) sites of TO ~ 0.742 /xg, in line with the experimen- 
tal value of m ~ 0.83 [Ib^^ 

Since the model is not frustrated, we expect that the 
above prediction from linear spin wave theory is a good 
approximation. To verify this conjecture, we simulate 
the static structure factor § corresponding to the prop- 
agation vector of the collinear state (7r,0) using QMC, 
which readily yields the order parameter S [Eq. (33) in 
Rcf. 58]: 

S\^) = (9) 

The results are presented in Fig. 10. We evaluated fi- 
nite 2D lattices with up to TV = 48x48 sites, and used 
Eq. (39b) from Ref. 58 for the finite-size scaling (Fig. 10): 

where 6*00 is the magnetic moment, extrapolated to the 
thermodynamic limit. The resulting 5*00 — 0.3297 yields 
'm = gfiBSoo — 0.717 ^B, corroborating the linear spin- 
wave theory result. 
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FIG. 10. (Color online) Finite-size extrapolation of the quan- 
tity S^{N) defined in Eq. (9). Finite- A'' values (crosses) were 
evaluated using QMC simulations for the effective model of 
the Cu(2) spins for a series of finite lattices up to A'^ = 48x48 
sites. The fit to Eq. (10) yields cri= 0.5509, a2 = 0.6436, 
CT3 = —0.5097, and Soo — 0.3297, with the extrapolated mag- 
netic moment m — gfisSoo — 0.717 fiB- 



C. Staggered polarization on the Cu(l) dimers 

The effective low-energy theory derived above in Sec. 
V A is a theory for the projection of the full wave function 
of the system onto the manifold where all Cu(l) dimers 
form singlets. However, there is also a finite GS compo- 
nent out of this manifold. In particular, once the Cu(2) 
sites order magnetically in the collinear stripe phase (see 
right panel of Fig. 8), they will exert a finite exchange 
field on the Cu(l) dimers. A simple inspection of Fig. 11 
tells us that this field is staggered, i.e., the local fields on 
the two Cu(l) spins of any given dimer arc opposite to 
each other. Specifically, the fields hi_2, as shown in Fig. 
11, arc given by 

hi,2-± (Jti - Jt2 + Jit)(Scu,) = T0.55Jd(Scu,), (11) 

where the last equation follows again from the DFT- 
based estimates for the ratios of the leading couplings 
(Table II, first row). 

Now, the important point is that the staggered field 
does not commute with the exchange interaction on the 
dimer, and as a result it can polarize the system imme- 
diately. This is very different from the case of a uniform 
field, where one must exceed the singlet-triplet gap Jd in 
order to (uniformly) polarize the sites of the dimer. 

To assess the amount of the induced polarization on 
the Cu(l) dimers, we consider a single AFM dimer with 
the exchange Jd in the presence of a staggered field h, 
described by the Hamiltonian 



H = Jd Si • S2 - h{St - SI) 



(12) 



The eigenstates of the dimer for h ~ are the singlet 
GS, |s) = (Iti) ^ lit)): s-nd the triplet excited states 

\h) = itt), \t-i) = III), and i^o) = ^(in) + i;t)). 

The staggered field leaves j^i) and |i-i) intact, namely 




OCu(1) 
• Cu(2) 



FIG. 11. (Color online) The classical ground state of the 
effective model derived in Sec. VA. The vectors hi, 2 denote 
the local exchange fields on the Cu(l) sites exerted by the 
magnetically ordered Cu(2) sites (blue and red arrows). 



T~l-\t±i) = -fl^ii), but it mixes |s) with \to) as follows 



n\to} ^ •■^\to} - h\s} 
n\s)^-^\s)^h\to) 



(13) 
(14) 



A straightforward diagonalization of the Hamiltonian in 
the subspace of {|s), |to)} leads to the eigenvalues 



(15) 



For h = 0, these results reduce to the energies Jd/4 and 
— 3Jd/4 as expected. The GS is given by 




1 

N 



+ \s)+h\t„) 



(16) 



where = 2h'^ + + Ja^JI/A + K^. As expected, 
|V'-)i->|s) for /ii-^O, while for /i > Jd, |^/'-)h> \\l). 

We may also look at the GS expectation value of the 
local polarizations on the dimer. We find 



2h^ + 4/2 + J^^Jl/ITl^ 



(17) 



For h = we recover jS'f 2I'/'-) = 0, while for h ^ Jd 
wc get maximum polarizations ("i/)- |S'f 2|'0-) il/2, 
which corresponds to the state | fi) ■ In the linear re- 
sponse regime h-^J^, we get (i/'_|5'f_2|V'-) — ^h/J^, 
thus the staggered susceptibility is Xs = V'^d- 

We can now apply these results to our case. Using 
Eqs. (11) and (17), the polarization on the Cu(l) sites 
can be estimated as 



K^cuJI ^0.176 ~0.516|(5cu,)| 



(18) 
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which is again in very good agreement with the experi- 
mental |(S'cui)|=0.54|(S'cu,)| from Ref. 26. 

D. The 1/2-plateau phase and the critical field Hci 

The separation of energy scales in CdCu2 (603)2 is a 
natural reason for the formation of the 1 /2 magnetization 
plateau. Due to the strong Jd, the Cu(l) spins involved in 
the Jd-dimcrs arc much harder to polarize with a uniform 
field than the weakly coupled Cu(2) spins. Naturally 
then, the 1/2-plateau corresponds to the state with the 
Cu(2) spins fully polarized and the Cu(l) dimers still 



carrying a strong singlet amplitude and a finite staggered 
polarization. 

Here, we go one step further and calculate iJci, the 
onset field of the 1/2-platcau by considering the one- 
magnon spectrum of the effective model in the fully po- 
larized phase. The one-magnon space consists of states 
with all but one spins pointing up. We label the possible 
positions of the down spin by one index for the unit cell 
and another index a = 1 — 4 which specifies one of the 
4 inequivalent Cu sites within the unit cell (sec sites 1-4 
in the right panel of Fig. 8). A straightforward Fourier 
transform gives the following 4x4 Hamiltonian matrix 
in the one-magnon space 



'^^l-magnon(k) — — 



2A 
J 



J 
2A 



Jyil 



J,(l + Jyie-"'^ 

Jy{l + e-'''y) J^{e-'^y 

2A 

J 



_|_ Q-i{k^+ky)-^ \ 

_j_ g-i(fcx+fcy)^ 

J 

2A / 



where we measure all energies relative to the energy of the 
fully polarized state, and A = gfJ-sH — (J^ + Jy + J/2). 
Using Eqs. (2)- (4), we find that one of the four branches 
of the one-magnon spectrum becomes soft (reaches zero 
energy) at k = when gfisH = J + 2Jx = ^^Jd- This 
soft mode signals the instability of the fully polarized 
state, hence 



H, 



cl 



1109 Jd 
6400 5^ 



(19) 



With 5~2.175 and Jd = 178K (AMF solution from Ta- 
ble II), we obtain iJci — 21.11 T, which is in very good 
agreement with the experimental value of Hd = 23.2T. 
This confirms that the third-order perturbation theory 
gives an adequate quantitative description of the low- 
energy physics. It also confirms that our set of exchange 
parameters gives an accurate quantitative description of 
the magnetism in this compound. 



E. Exact Diagonalizations 

We are now going to check the above physical picture 
obtained from the effective model by numerical exact di- 
agonalizations on the original model. We first address 
the nature of the magnetic GS by evaluating the spin- 
spin correlations (Sq-Sr) as a function of the "distance" 
R (Fig. 12). In order to elucidate the difference in the 
coupling regimes of Cu(l) and Cu(2) spins, we choose 
two independent paths on the 2D spin lattice, with dif- 
ferent initial spins So (Fig. 12, lower panel). As ex- 
pected for strong Jd , the structural dimers exhibit strong 
AFM correlation (Sq-Si) = -0.57 (green line in Fig. 12), 
while the Jit bonds bear much weaker FM correlation 
(So • Si) =0.15 (brown line in Fig. 12). 



A rough estimate for the ordered magnetic moment 
can be obtained from the spin correlations of the maxi- 
mally separated spins on a finite lattice. The small size 
of the 32-site finite lattice impedes a direct evaluation of 
the order parameters. However, the same-sublattice cor- 
relations for the Cu(l) and Cu(2) sublattices, evaluated 
using the same N = 32 cluster, can be compared to each 
other. In this way, finite size effects, which are crucial 
for the values of tocu(i) and 'm,cu(2), should be largely 
remedied for the mcu(i):TOcu(2) ratio. 

The maximal separation in terms of exchange bonds 
between the same-sublattice spins on the 32-site fi- 
nite lattice is five and corresponds to (So ■ Ss) corre- 
lations in Fig. 12. ED yields (Sq • Sg) [Cu(l)] = 0.060 
and (So • S8)[Cu(2)]= 0.174 for Cu(l) and Cu(2), 
respectively. In the simplest picture, the quan- 
tity ^(So-S8)[Cu(l)] : ^(So-S8)[Cu(2)] =0.59 should 
be close to the ratio of the magnetic moments. Al- 
though the results may be still affected by finite- 
size effects, a comparison to the experimental^^ 
™Cu(i) '■ "T'Cu(2) — 0.45 : 0.83 = 0.54 reveals surprisingly 
good agreement between theory and experiment. 

Next, we examine the evolution of the GS expecta- 
tion values of the local spin-spin correlations in magnetic 
field, which are presented in Fig. 13. For h* = 0, the cor- 
relations within the structural dimers and along the Jt2 
bonds are dominant, amounting to —0.570 and —0.366, 
respectively. The first one is somewhat weaker than the 
pure singlet value of —0.75, because of the finite admix- 
ture with the triplet \to) in GS, as explained in Sec. VC. 

Now, according to Fig. 13, a magnetic field enhances 
the intradimer correlation. This counter-intuitive result 
can be physically understood as follows. In a finite field, 
the Cu(2) spins order in the plane perpendicular to the 
field (this happens in an infinitesimal field, since the crys- 
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FIG. 12. (Color online) Upper panel: spin correlations 
(So ■ Sr) in the magnetic ground state of CdCu2(B03)2, sim- 
ulated on the 32-site finite lattice with periodic boundary con- 
ditions using ED. Lower panel: A'' = 32 sites (full plot) and 
A'^ = 16 sites (shaded rectangle) finite lattices used in our ED. 
Sq-Sr pathways are used for the (Sq-Sr) (R) plot in the upper 
panel. For the notation of exchange couplings, see Fig. 5. 



talline anisotropy effects are neglected) and at the same 
time develop a uniform component along the field. While 
the former, as discussed in Sec. VC, gives rise to a stag- 
gered exchange field on the Cu(l) dimers, the latter in- 
duces a uniform exchange field. By increasing the ex- 
ternal field, the staggered component decreases, while 
the uniform component increases. At the 1/2-plateau, 
the Cu(2) spins are almost fully polarized, hence only 
the uniform component survives. Now, in contrast to 
the staggered component, the uniform field commutes 
with the exchange interaction on the Cu(l) dimer and 
thus leaves the singlet GS wave function intact. In other 
words, the magnetic field suppresses the admixture of a 
triplet in the singlet GS. The steep decrease in the cor- 
relation along the Jt2 bonds in a field reflects the same 
physics. 

At the critical field /i*~0.2 (24 T), GdCu2(B03)2 
enters the 1/2-plateau phase. Here, in contrast with 
the strong AFM correlations within the Cu(l) dimers 
((So • Si) = —0.718), the correlations along the Jti, Jt2 
and Jit bonds are much smaller and similar in size, but 




h* 

FIG. 13. (Color online) Magnetization (5?) for Cu(l) and 
Cu(2) spins, as well as spin correlations (S; ■ Sj) correspond- 
ing to the leading couplings Jd, Jti, Jt2, and Jit (line styles 
according to the inset) as a function of magnetic field h* in 
units of Jd. 



have different signs. The negative correlation for Jt2 and 
the positive for Jit are in accord with the nature of these 
two couplings, AFM and FM, respectively. For weaker 
AFM Jti , the correlation is positive due to the frustrated 
nature of the spin model. 

Since the local spin correlation on the Cu(l) dimers is 
very close to the singlet value, the contribution to the 
total magnetization from the Cu(l) spins is expected to 
be negligible. Therefore, the magnetization results pre- 
dominantly from the Cu(2) spins. This picture is sup- 
ported by simulations of the local magnetizations (5*1), 
shown in the upper panel of Fig. 13. In the 1/2-plateau 
phase, the local magnetization of Cu(2) is (S'^^j) — 0-483, 
while the magnetization of Cu(l) is negligibly small with 
i^Cni) =0.017, which also agrees with the main picture 
obtained from the effective model. Therefore, within the 
error from finite-size effects, only about 3 % of the Cu(l) 
moments are polarized in the 1/2-plateau phase. 

Let us now discuss the width of the 1 /2-plateau. This 
corresponds to the energy we need to pay to fully polar- 
ize a Cu(l) dimer, i.e. to convert a singlet into a triplet 
excitation. The width of the plateau can be estimated by 
neglecting the pcrturbative corrections that give rise to a 
hopping of these triplet excitations (and, thus, to a small 
dispersion). Our ED results (upper panel of Fig. 13) 
agree very well with this estimate: the critical field at 
which the 1/2-plateau ends is h* = 1.17 (143T). 

Now, once we pay the energy Jj to excite a triplet, 
we can also polarize the remaining Cu(l) dimers very 
quickly due to the small (pcrturbative) kinetic energy 
scale of the triplet excitations. Indeed, our ED results 
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for the magnetization above Hc2 show a rapid increase 
up to the saturation value. In addition, we find no indi- 
cation for any other magnetization plateaux above 1/2, 
which suggests that the interaction between the triplet 
excitations is negligible. 



VI. SUMMARY AND OUTLOOK 

The magnetism of CdCu2 (603)2 was investigated on 
a microscopic level by using extensive DFT band struc- 
ture calculations, ED and QMC simulations, third-order 
perturbation theory, as well as high-field magnetization 
measurements. In contrast to the previously suggested 
"chains + dimers" model, we evaluate a quasi-2D mag- 
netic model with mostly antiferromagnetic couplings: an 
intradimer exchange Jd and two intra-tetramcr interac- 
tions Jti and Jt2, as well as a ferromagnetic coupling 
Jit between the tetramers. The Jd : Jti : Jt2 ■ Jit ratio is 
close to 1 : 0.20 : 0.45 : —0.30 with a dominant Jd of about 
178 K. Topologically, this microscopic model can be re- 
garded as decorated anisotropic Shastry-Sutherland lat- 
tice (Fig. 1). 

While the spin pairs forming the structural dimers 
show a clear tendency to form magnetic singlets on the 
Jd bonds, the low-temperature magnetic properties of 
CdCu2 (603)2 are governed by interdimer couplings. In 
contrast to the strongly frustrated SrCu2(B03)2 featur- 
ing a singlet ground state, the inequivalence of Jti and Jt2 
largely lifts the frustration and triggers long-range anti- 
ferromagnetic ordering. This model correctly reproduces 
the substantial difference between Cu(l) and Cu(2) or- 
dered moments, as observed experimentally. Our simula- 
tions and the subsequent fitting of the magnetic suscepti- 
bility and magnetization further support the microscopic 
model and the proposed parameterization. 

It is interesting to compare the key results for 
CdCu2 (603)2 with the archetype Shastry-Sutherland 
S = ^ system SrCu2(B03)2. Although the substitution 
of Sr by Cd significantly alters the crystal structure and 
stabilizes a different magnetic ground state, the micro- 
scopic magnetic models of both systems bear apparent 
similarities (Fig. 1). Moreover, the dominance of the Jd 
coupling in CdCu2(B03)2 leads to an effective nonfrus- 
trated model for the Cu(2) spins, which is topologically 
equivalent to an anisotropic Shastry-Sutherland lattice 
model. It is noteworthy that the phase diagram of such 
a model was recently discussed. 

One of the main objectives of this study is to stimulate 
further experimental activity on CdCu2 (603)2. Such 
studies are necessary to refine the parameters of the sug- 
gested microscopic model and address yet unresolved is- 
sues. Our extensive analysis of the magnetic properties of 
this complex system delivers several results that should 
be challenged experimentally. For example, we have re- 
solved that the Cu(l) dimers develop a sizable staggered 
polarization which is suppressed by an applied field and 
vanishes at the 1/2 plateau. This behavior could be ad- 



dressed experimentally by nuclear magnetic resonance. ^'^ 
On the other hand, inelastic neutron scattering experi- 
ments could challenge the magnon dispersions calculated 
for the effective low-energy model. In addition, the mag- 
netic couplings, associated with edge-sharing and corner- 
sharing connections of CUO4 plaquettes, are strongly de- 
pendent on Cu-O-Cu angles. Therefore, the magnetic 
properties of CdCu2 (603)2 could be tuned by applying 
external pressure. Similar experiments on SrCu2 (603)2 
(Ref. 11) reveal excellent potential of such studies. 
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Appendix A: Linear spin wave theory around the 
stripe phase of the effective model 

Here, we provide the details of the linear spin wave the- 
ory around the stripe phase of the Cu(2) spins. We begin 
with the effective model Hamiltonian that was derived in 
Sec. VA for the Cu(2) spins: 

'H = J (^Si,u • 82,11 + 83,11 • 84,11^ 

u 

+ Jx 83,11 • (^Si u + 8i,„4.t^^ 

+ Jy 83,,, • ^82, u + 82, u+t„^ 

Jx + Jx 82 

+ Jy 81, „ • S4^u-t^~ty + Jy 84, „ • 8i,„-)_j^ (A-1) 

where tx,y = 2a e^^y are the primitive translation vectors 
in the model (a is the lattice constant of the square lat- 
tice), and the index u labels the unit cell. The remaining 
spin operator index (1-4) labels the 4 sites of the unit cell 
(see right panel of Fig. 8). We then take as our reference 
classical state the one where the spins 1 and 3 (resp. 2 
and 4) of the unit cell point along the positive (resp. neg- 
ative) z-axis, and perform a standard Holstein-Primakoff 
expansion in terms of four bosonic operators labeled as 
c„ , dy^ , which correspond to the four spins per unit 
cell S'i,u, S'2,u, 5'3,u, Si^u respectively. The transformation 
reads 



'-'l,u 


~ -s 


+ at an, 


'-'l,M 


^ V25a+ 


QZ 


~s- 




Q+ ^ 
■^2,11 — 


^25 bu 


QZ 


~s- 




Q+ ^ 
'-'3,11 — 


V2Scu 


qZ 
'-'4,11 
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"-'4,'U 


^V2Sdt 
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Replacing these expressions in the Hamiltonian and per- pact matrix form 
forming a Fourier transform, we get 



S 



n = Eob+H2b + 0{S'') 



(A3) 



n2b^6E, + -J2K ■'^k-Ak (A4) 

k 

where SEi = -~2^NucS, 

Afc=(a^, b-^, 4, dl, a-k, b-k, C-k, d-fc ) 



where Eob = —"^^^.N^cS^ stands for the classical energy, 
and ^ = J + 2Jx — 2Jy. The term 7^26 stands for the 
quadratic Hamiltonian which can be written in the com- and 



Mi 
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Jye-''''^{l + e-''=y)\ 

Jyil + e-^'y) 




\Jye'^'^{l + e'''y) 
-J 

-J 

'J^l + e'''-) 

-7^6'*=" (1-1- 6^*=-) 
















-J 
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(A5) 
(A6) 



(A7) 



(A8) 



To diagonalize 7^26 we search for a new set of bosonic 
operators Afc given by the generalized Bogoliubov trans- 
formation Afe = Sfc- Afc, such that the matrix S^M^Sfe = 
Jlfe becomes diagonal. The transformation must also pre- 
serve the bosonic commutation relations, which can be 
expressed compactly as g = g = Sa_, • g • S^, where g is 
the "commutator" matrix 

g - A,. . A+ - ((A^)" • A^)" = ( -1, ) (A9) 

and I4 stands for the 4x4 identity matrix. The above 
two conditions give 

(gMfe) • Sfc = Sk ■ (gflfc) = Sfe • n'k (AlO) 

which is an eigenvalue equation in matrix form (the 
columns of Sfc contain the eigenvectors of gMfc). 

One can show^^ that if Mfc is semi-definite positive, 
then the eigenvalues of gMj, come in opposite pairs: 

where uifc is a diagonal matrix with non-negative entries 
Wife,. . ., aj4fe. This in turn leads to 

'H2& = 6E1 + 6E2 + 5*^ (ujik a^dk +■ ■ .+ W4fe d't'^kj , 

(A12) 



where the second correction 

dE2 = ^J2 ^^^^ +■ ■ •+ ^4fc) , (A13) 

k 

stands for the total contribution from the zero-point en- 
ergy of all independent harmonic operators involved in 
the theory. 

1. Renormalization of the magnetic moments 

Let us now look at the effect of the quadratic fluctua- 
tions on the magnetic moments. We consider the spin Si 
operator inside the unit cell n = 0. The classical vector 
points along the z-axis. We have 

k,q 

Using the transformation A/^ = • we get 

(A15) 

In the vacuum GS, the only non- vanishing expectation 
values are of the type {d^d'^) = 1. Thus from the above 
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sum we should keep only terms with z = j = 5 — 8 and k = q, namely 

(SD=5-^5]^|5.(l,z)r (A16) 

k 4=5 

The second term can be calculated by a numerical inte- 
gration over the BZ. 
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